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Abstract. Accretion onto a compact object must occur through a disc 
when the material has some initial angular momentum. Thin discs 
and the thicker low radiative efficiency accretion flows (LRAFs) are 
solutions to this problem that have been widely studied and applied. 
This is an introduction to these accretion flows within the context of 
X-ray binaries and cataclysmic variables. 

00 

Accretion describes the increase of the mass of an astrophysical object when 
fNJ i matter is deposited onto it. The conversion of a fraction of the gravitational poten- 
tial energy of the accreted material into radiation can dominate the emission from 
i the system. Accretion is a major source of energy in a varied lot of astrophysical 
objects including interacting close binaries, active galactic nuclei and protostcllar 
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systems (the standard textbook is Frank et al. 2002). 

My focus is on close binaries in which a normal star transfers mass onto a black 
hole/neutron star (X-ray binaries, XRBs) or white dwarf (cataclysmic variables, 
O | CVs) via Roche-lobe overflow. Exhaustive reviews of their observational proper- 



ties may be found in Lewin et al. (1995| ) and Warner (1995) respectively. In these 
systems, an accretion disc forms around the compact object due to angular mo- 
mentum conservation. The detailed characteristics of the flow depend upon how 
matter dissipates its potential energy and gets rid of its angular momentum in the 
disc. This basically leads to two extremes: the "low radiative efficiency" accretion 
?H ' flows (LRAFs, which are geometrically thick) and the radiatively efficient, geo- 
metrically thin discs (hence the title). This is an introduction to these solutions. 
§1 groups the basic tools needed to present the solutions in §2 (thin discs) and §3 
(thick discs); §4 is a brief discussion of the transition from one to the other. 



c/5 



1 The building blocks 

1.1 Accretion efficiency, Eddington luminosity 

In Newtonian dynamics, the available gravitational energy from a small mass m 
moved from infinity to the surface i?* of the accretor is GM+m/ R+. The energy 
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released can be a sizeable fraction 77 of the rest-mass energy mc 2 when the accretor 
is a neutron star or black hole: assuming a thin disc and free fall at the inner- 
most stable circular orbit then r\ w 0.42 for a maximally rotating Kerr black hole 
(Novikov & Thorne 1973; see also Gammie 1999). 

The Eddington luminosity Le is the luminosity for which radiation pressure 
exactly balances the gravitational pull on the accreted material. Accretion is 
quenched above this luminosity (but see Begelman 2001; Shaviv 2001). For a 
spherical inflow of ionised hydrogen such that the radiation pressure is due to 
Thomson scattering on electrons and the dominant gravitational pull is on the 
protons, Le is: 



R 2 



Le ot 
4itR 2 ~ 



hence L E » lO 38 (M*/M©) erg • s" 



(1.1) 



The luminosities of XRBs and CVs are generally lower than this limit ( Lewin et al 



1995; Warner 1995 ). The Eddington mass accretion rate Me is defined as: 



# E c 2 = L E i.e. M E « 10 18 (M*/M )(O.l/r?) g s" 



(1.2) 



1.2 Accretion flow temperature 

Assuming the gravitational energy is released into radiation at the surface of the 
compact object, the minimum temperature of the flow is that of the black body 
radiating the same luminosity. At the Eddington limit this gives: 

47ri? 2 crT4 = L E i.e. kT bb « 1.5 keV (L E /10 38 erg s -1 ) 1/4 (.R*/10 km)" 1 / 2 (1.3) 

Accreting black holes and neutron stars should radiate mainly in soft X-rays while 
white d warfs (R+ ss 10 4 km) should radi ate in UV. This is consistent with obser- 
vations ( Lewin et al. 1995 ; Warner 1995 ). 

On the other hand, the maximum temperature is obtained when all the grav- 
itational potential energy is transformed into thermal energy e without radiation 
losses (adiabatic flow): 

e = 5_ = — — i.e. kT g « 45 MeV (MJM Q )(R*/10 km)" 1 (1.4) 

2 /imn it* 

assuming ionised hydrogen (/i = 0.5). This is twice the virial temperature of bound 
particles in circular orbit at i?*. Accretion onto a compact object can power the 
emission of gamma rays. 



1.3 Disc formation 

Matter infalling onto the compact object will generally have some non-zero angu- 
lar momentum. A particle with a ballistic trajectory will therefore have R 2 Q = 
(i? 2 f2)t=o- At the radius of closest approach to the compact object l/2(i?£!) 2 = 
GM/R C so that 2GMR C = (R 2 ft)l. Only a particle with very low angular mo- 
mentum can directly hit the compact object (R c < i?*). This condition is not 
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met in compact binaries where matter comes from an orbiting companion and the 
stream will go round the compact star and intersect itself. Subsequent shocks lead 
to the dispersion of energy and the stream settles onto a circular orbit with the 
initial angular momentum (see Lubow & Shu 1975 for details). 

A steady supply of matter with some specific initial angular momentum piles 
up in a ring at this circularisation radius. There is no accretion unless some matter 
can transfer its angular momentum to reach smaller orbits. Under such a process, 
an accretion disc forms extending down to the compact object. 



1.4 Angular momentum transport 

One process by which particles may exchange angular momentum is viscosity (see 



e.g. Terquem 2001 for an introduction). In gas kinetic theory, molecular viscosity 
arises from the exchange of momentum across the surface of a fluid element. The 
resulting force is proportional to ^ mo i ~ Xu where A is the mean free path between 
collisions and u is the mean thermal speed of the particles. Assuming particles in 
a plasma interacting only via Coulomb forces leads to: 

A ~ 1 cm (T/10 5 K) Go/10" 8 g cm" 3 ) 
u ~ 10 6 cm s" 1 (T/10 5 K) 1/2 



using typical accretion disc temperatures and densities ( Frank et al. 2002 ). The 
viscosity coefficient is therefore i/ mo i ~ 10 6 cm 2 s -1 . 

Anticipating a little bit on result from the thin disc model, we can get an 
order-of-magnitude observational estimate for visc osity in accretion discs. If the 
eruptions of dwarf novae, a sub-class of CVs (see j j2.8| ), are due to the accretion 
of matter coming from the outer regions of a disc and transported by viscosity 
then t'disc ~ R<iv r where i?d is the outer disc radius and v T is the radial velocity 
of the infalling matter. v r ~ Rd/Te where r e is the timescale of the eruption. 
The disc radius can be estimated using e.g. the circularisation radius (§1-3) and 
typically will be ~ 10 10 cm. For an eruption lasting a day, ^di sc ~ 10 15 cm 2 s _1 
i.e. fdisc 3> Vmo\- Molecular viscosity is much too weak to account for accretion in 
discs. 

The process by which angular momentum is transported has been the subject 
of intense research with turbulent transport the prime suspect. Turbulent viscosity 
is often modelled using the Navier-Stokes formalism but with a coefficient i^urb ~ 
Aturb^turb where the scale and speed are those of the turbulent eddies. These are 
dynamic properties of the fluid, not intrinsic as with molecular viscosity. In a disc, 
the size and speed of the eddies are likely to be limited by the scale-height H and 
sound speed c s : 

v = ac s H (1-5) 
where a is a parameter < 1 . This is the famous a parameterisation of disc viscosity 



first described in a landmark paper by Shakura & Sunyaev (1973). In the following 
I will use this parameterisation. 
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One would rather want to derive v from first principles. At present, only tur- 
bulence arising from the magneto-rotational instability (MRI) has succeeded in 
predicting any significant viscosity (see Balbus & Hawley 1998 for a review). For- 
tunately, there is some sense in modelling the angular momentum transport and 
energy dissipation of MHD turbulence using the a parameterisation ( Balbus & Pa 
paloizcu 1999). However, the mechanism may not work in weakly ionised flows and 



this is a problem for models with a cold accretion disc (Stepinski et al. 1993; Gam- 
mie & Vlenou 1998). Putative global hydrodynamical instabilities might take over 
in these conditions. There are other possibilities for angular momentum transport 
including instabilities in a self-gravitating disc (e.g. Balbus & Papaloizou 1999 and 



references therein), spiral waves excited by tidal torques (Spruit 1987), hydromag- 
netic winds launched from the disc ( Blandford fc Payne 1982) ) , radiative viscosity 
(e.g. Loeb & Laor 1992) etc. However, these only apply in specific conditions. 



1.5 Vertically integrated disc equations 

The disc equations derive from the equations of fluid dynamics combined with a 
model for viscosity (which contains the turbulent magnetic field contributions in 
the MRI case), an equation of state for the gas and a description of the radiative 
processes. A standard set of assumptions to start with is: 

• Axisymetry so that d/d<j) = in cylindrical coordinates. 

• The only non-zero component of the stress tensor is the azimuthal shear r r ^. 

• The gas is perfect so the internal energy per unit mass is e = c v T and the 
gas pressure P = pe{^ — 1) with 7 = c p /c v . Radiation pressure is neglected 
here for simplicity. 

In addition I suppose the disc self-gravity is negligible (this assumption must be 
dropped for protostellar discs and AGNs) and that relativistic corrections are neg- 
ligible (which is fine when more than a few gravitational radii away from the 
compact object). Assuming hydrostatic balance so that v z = 0, the vertical mo- 
mentum conservation is [P = pc 2 s ): 



dP 

Oz 



-P9z 



d\nP 
d In z 



fl 2 



1 



z 

R 2 



(1.6) 



If the height H from the midplane of the disc is < R high order terms in z/R 
can be neglected. This equation can be integrated analytically for a perfect gas 
yielding P(z) and p(z). Averaging P over z gives a relationship between H and 
the mid-plane sound speed which, to a factor of order unity, is : 



H = Cs/Qi 



(1.7) 



The scale- height c s /Ok appears in Eq. 1.6. The detailed vertical balance can be 
bypassed by assuming the disc height is given by this relationship. The Shakura- 
Sunyaev prescription for the viscosity becomes 



ac 2 jn K 



(1.8) 
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Note that v is effectively integrated over z. With a Navier-Stokes formulation of 
the viscosity, the integrated stress is: 

on 



dR 



where E = 2p H is the column density and p a is the mean density. The radial 
evolution equations are then integrated and decoupled from z, resulting in a set 
of time-dependent ID equations: 



9E 
~dt 



R~dR 

dv r 



dR 2 n 



at 

ds 

T di + VrT dR 



(ERv r ) 
dv r 

- V r 

dR 
dR 2 fl 
r OR 
ds 







Rn 2 - RCl^ 



1 dP 
Po dR 



1 d 
Ei? <9i? 
1 



r M 2 r, 



r&) 



1 d 
Ei? 9i? 



dR 



m ls -Q 



rad/ 



(1.9) 
(1.10) 

(1.11) 
(1.12) 



In order, we have the mass, radial momentum, angular momentum, and energy 
conservation equations. Rfl and v r are the angular and radial velocities while 
P = p c 2 - c s and therefore T are the mid- plane values of the sound speed and 
temperature in this framework (Eq. 1.8). Again, only the dominant term in the 
radial component of gravity is kept g r ~ RQ^. In the thermal equation, s is 
the entropy and Q^ is is the energy generated by the work of viscous forces in the 
annulus: 



^E R 



Qmd represents radiative losses. Using ds 
can be rewritten as: 



OR 

de 



(1.13) 

Pd(l/p ) the thermal equation 



de 



de 



dt +Vr 'dR 



P 



dp 
dt 



dp 
dR 



q; 



Ml J 



and e can be written using T g (see Eq. 1.4): 

,2 T 



(7-1) 



(1.14) 



(1.15) 



Other versions of these equations accommodate additional energy sources and 
sinks, conduction or convection terms, radiation pressure, different proton and 
electron temperatures, magnetic fields, relativistic corrections etc. 



1.6 Steady state disc 

A disc in stei 
over R gives: 



A disc in steady state has Jt = 0. Combining Eqs. 1.9 and 1.11 and integrating 



R 2 Q 



vR 2 dfl 
v r dR 



constant 



(1.16) 
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Assuming there are no viscous torques at the surface of the compact object, the 
integration constant is R^Cl*. This term is very small compared to R 2 Q, far from 
the boundary layer. Writing 6=1 — (i?^fi*/i? 2 fi), Eq. 1.16 then gives 



bv r 



v dhn.fl 
Rd\nR 



(1.17) 



The radial momentum equation can be rewritten using Eqs. 1.17 and 1.15 as: 



9l. 

n 2 



1 = ( 7 - 1) 



T\ dlnP 



T e J dlnR 



+ (7 -I) 5 



a T 



d In £1 \ d In 



dlnR J dlnR 



(1.18) 



Defining the cooling efficiency / as 1 — Q ra( j/Qvis> the energy equation becomes: 



d\nR 



d In e . . d In p Q 

-(7-1)' 



dhxR 



dlnR 



(1.19) 



In a solution to the equations, / should be calculated self-consistently from the ra- 
diative processes. I find this set of dimensionless equations useful when discussing 
the assumptions that go into a thin or thick disc. 



2 Radiatively efficient flows: thin discs 



Thin accretion discs are flows in which the energy liberated by accretion is effi- 
ciently radiated locally. The next section (§3) deals with the LRAFs whe re this is 
not the case. Thin discs were first described by Pringlc fc Rees (1972 ) and ghakura 
fc Sunyaev (1973). The standard presentation of t his model may be found in the 
review by Pringlc (1981 , see also Frank et al. 2002 ). 



2.1 Main properties in steady state 

For a cool steady disc such that T/T g <C 1 then we must have / ~ i.e. Q^ is w 
<3Z.j, viscous heating is balanced locally by radiative losses. In other words, the 
disc is radiatively efficient. The definition Eq. 1.15| of e as a function of T/T g 
shows that c s i?f2 K hence 

H/R<^1 

A cool accretion flow is geometrically thin. It is easy to see from the dimensionless 
equations that T/T g -C 1 further implies that a steady flow is Keplerian with 
f2 ps fix, azimuthally supersonic with ~ RQk 3> c s but that the radial inflow 
velocity 

6«r = -~<C a (2.1) 
z K 

is highly subsonic. The mass accretion rate in the disc, defined as 



M = -inRT.Vr 



(2.2) 
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is a constant (Eq. |LS| ). The balance of heating and cooling gives: 

Q+ = -vEfi* =Q- = 2aT c 4 ff (2.3) 
The factor 2 accounts for radiation losses from both sides of the disc. This can be 



rewritten using Eqs. 2. 1 - 2.S as 



4 3b GM+M 

aT ^ = 8^^r- ( 2 - 4 ) 

The disc dissipates an energy 3/2GMM/R 2 dR between R and R + dR far from 
the inner boundary, three times more than the local release of gravitational energy. 
The extra energy comes from smaller radii with viscous transport redistributing the 
way the available gravitational energy is radiated in the thin disc. This equation 
also shows clearly that the emission of a steady state disc is independent of the 
particular heating mechanism (viscosity in this case). Thus, there is no need to 
know the local angular momentum transport mechanism in order to predict the 
disc luminosity. 

Eclipse mapping has been used to test the radial dependence of the emission. 
In a high inclination system, the companion gradually eclipses different parts of 
the disc during its orbital motion. The resulting eclipse profile can be used to 
map the radial distribution of the disc brightness temperature. Results show that 
dwarf novae in outburst and persistent CVs (novae-like systems) have temperatures 
varying as i? -3 / 4 in agreement with expectations. On the other hand, quiescent 
dwarf novae show flat radial temperature distributions and that supports the disc 



instability model which is discussed later in § 2.7 (not all observations fit this simple 



picture, see e.g. 3mak 1994; Warner 1995; Baptista 2001) 



2.2 Vertical structure 

The thermal balance (Eq. |2.3| ) is written without any knowledge of the radiative 
processes that cool the disc. This is needed to relate Q~ d to the quantities appear- 
ing in the radial equations. For instance, one can solve for Q~ by integrating the 
hydrostatic balance (Eq. |L~6| ) and vertical radiative transfer equation (Eq. p.. 12 ) 



explicitly. For optically thick radiation, this is like solving a stellar atmosphere (see 
Hubeny 199C| , and references therein). The caveat is that additional assumptions 



have to be made on how viscous dissipation occurs locally in the layer (typically 

The simple one zone model for the vertical structure has HQ^ = c s , S = 
2poH, t — kT, and P = pokT/pmn- The radiation pressure is neglected and k 
is the Kramers opacity (Rosseland mean of bb, bf & ff processes). If the layer is 
optically thick, i.e. the disc cools at the maximum possible rate, the Eddington 
approximation for the radiative flux relates T c g to T . These assumptions can 
be shown retrospectively to be valid for the range of densities and temperatures 
considered. 
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This completes the set of equations and a steady solution can now be found (Frank 



et al. 2002): 



S = 5.2 a-WM^M}'*^*'* g cm" 

H/R = 0.02 a-^M^M^R^ 8 

T = 14000 a" 1 / 5 MU W Mi m Ri m K 

t = 190 a-^ 5 MU 5 

v r = 3 • 10 4 a^M% W M^ /A R^ cm s 



2 — 1 

cm s 



where M\ is the mass of the compact object in solar units, Rio is the radius in 
units of 10 10 cm and M\q is the accretion rate in units of 10 16 g s" 1 . Around 
a black hole or neutron star the disc can extend down to small radii. As the 
temperature increases, electron scattering takes over from free-free as the major 
source of opacity. At very small radii and high accretion rates, radiation pressure 



can become dominant but the disc is then unstable (see §2.6) 



2.3 Radiation spectrum 

The outgoing spectrum of the optically-thick model is calculated by considering 
each disc annulus radiates as a blackbody at a temperature T g and integrating 



over R. The resulting spectrum is a modified disc blackbody (Fig. 20 of Frank et al 



2002 ) . The maximum temperature of the radi ation is reached close to the surface 
of the accreting object at R = (7/6) 2 i?* (Eq. fyj). oT c 4 ff is of the order of 1 eV 
in the case of white dwarfs, 10 eV for massive black holes and 1 keV for neutron 
stars and stellar mass black holes. The big blue bump in AGNs and the soft X-ray 
emission in X-ray binaries can be interpreted as modified blackbodies from a thin 
disc. Optical observations leave little doubt that thin discs are present in compact 
binaries (e.g. double-peaked emission lines, radial temperature profiles, statistics 
of eclipsing systems, orbital lightcurves etc). 

But it was quickly recognized that these objects also sh ow emission at higher 



energ ies, which cannot be explained by a thin disc (e.g. Lightman fc Shapiro 
1975| ). One possible origin for this radiation is the boundary layer where the an- 



gular velocity of matter decreases rapidly to match that of the compact object. 
Integrating Eq. 2.4 over R shows only half of the available energy GM+M / R% is 



radiated within the thin disc and the balance is dissipated in the boundary layer 
(when the object has a surface). Modelling the contribution of the boundary layer 



to the spectrum of CVs and neutron star XRBs is a difficult problem (see Regev 



1991; Popham & Sunyaev 2001, and references therein). Other likely possibili- 
ties include contributions from a tenuous hot corona above the disc (somewhat 
analogous to the solar corona; see Poutanen 199§| and references therein) and/or 



emission from a different type of accretion flow closer to the compact object (see 



2.6 and 3.2) and/or a jet (e.g. Markoff et al. 2001) 
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2.4 Evolving thin disc: timescales 

A cool evolving disc is not necessarily Keplerian and this assumption must be 
made when studying time-depend ent thin discs. Numeric al simulations show this 



is a very good approximation (e.g. Ludwig fc Meyer 1998 ). The radial momentum 



equation can be put in a dimensionless form: 

Rd\nv r din \ O 2 . . fT\dlnP n .„ „, 

R^\V r — +m^R) = tii- 1) \tJ 8^r kQ (2 - 6) 

showing v r <C RQk, just like the steady case. The conservation of angular mo- 



mentum with £1 = f^K leads to v r ~ v/R (see Eq. 1.17 ). 

The major timescales of a thin disc are directly identified from the dimension- 
less forms of the other Eqs. L^-1.12. The dynamical timescale is: 



idyn ~ — ~ 7^- « 100 s M x 1/2 R% 2 (2.7) 



This is the shortest timescale over which azimuthal inhomogeneities (e.g. flares) 
would be expected to vary. Assuming the typical scale for the radial variation is 
R, the angular momentum equation yields the viscous / accretion timescale: 

t vis ~ V ~ ( f ) td ^ * 4 dayS ^mTm^R 5 ^ (2.8) 

This is also the character istic time R/v r needed to accrete viscously matter ini- 
tially at a radius R (see § |2.5| ). This accretion timescale is much longer than the 
dynamical timescale. Perturbations to the hydrostatic balance propagate at the 
sound speed so the vertical timescale is: 

TJ 

tvert ~ ~ ^dyn (^'^) 

The characteristic time needed to establish hydrostatic balance is of the same order 
as the dynamical timescale i.e. short compared to the accretion timescale. The 
energy equation yields the thermal timescale: 

e c 2 1 

ither ~ —T ~ r^-? ~ — ^dyn (2-10) 

q+ uil K a 

This is the characteristic timescale needed to establish the local balance of heating 
and cooling and is short compared to the accretion timescale. 



2.5 Diffusion equation 



The mass and angular momentum conservation equations of a thin disc, rewritten 
in terms of S, can be combined into a single evolution equation: 
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This diffusion equation describes how t he column density E evolves from a given 
initial radial profile and is discussed in Lynden-Bell & Pringle (1974). 

The most pedagogical example is that of an infinitely thin annulus of matter 
which spreads both in and out under the influence of a constant viscosity v (Figs. 4- 
5 of Lynden-Bcll & Pringle 1974). After an initial adjustment, the mass accretion 
rate onto the compact object approaches a self-similar regime where M varies as 
some power of time t. This is characteristic of constant total angular momentum 
solutions. After a few viscous timescales t V i S ~ R 2 /v, most of the initial mass is 
accreted onto the compact object while a decreasing fraction of the matter carries 
away the angular momentum at infinity. In a binary, the thin disc has an outer 
radius set by tidal torques and angular momentum is given back to the companion 
(Papaloizou & Pringle 1977). Semi-analytical solutions to this diffusion equation 



under different assumptions can also be found in e.g. Bath & Pringle (1981) and 
Lyubarskii & Shakura (1987). 



2.6 Stability of steady thin discs 



A steady thin disc is suscep tible to thermal and viscous instabilities ( Shakura & 
Sunyaev 197(| |Pringlc 197^ ). The first arises when cooling cannot keep up with 



heating i.e. if heating increases faster than cooling when the temperature of an 
annulus is raised (the thermal balance is local). The stability criterion is: 



dQ+ 
dT 



fdQ- 



v < \9T 



(2.12) 



The thermal timescale is much shorter than the viscous timescale on which 
evolves (jjpH). Setting E constant, Q + ~ vE ~ T (Eq. Of). For optically thick 



radiative cooling then (Eq. 2.5) Q~ ~ T 4 /r 
opacity r ~ /9 T~ 7 / 2 E), hence the disc of 
medium is optically thin then Q 



2.2 



T i5/2 H /jV - T 8 ( ta k in g Kramers 
is thermally stable. However, if the 
1 (again for free-free radiation) 



and the annulus is unstable ( Pringle et al. 1975| ). 

The disc is viscously unstable when perturbations of the density are amplified 
rather than smoothed out by accretion. The criterion for viscous stability is: 



dM dvZ 
9E >0 ° r ^E >0 



(2.13) 



using M ~ vll (steady-state thin disc) . The criterion is easily found by linearizing 
Eq. 2.11 for small perturbations. Density perturbations grow on the long viscous 
timescale so thermal equilibrium can be assumed. Since Q + ~ i^Ef^ ~ PHUk 
and Q~ ~ T 4 /kE ~ P/E one has H oc 1/E in a radiation pressure / electron 
scattering (re constant) dominated annulus. Writing the hydrostatic equilibrium 
shows P ~ (E/lf)Cg ~ E77il^ oc 1 is independent of E hence vY, oc PH oc 1/E: 
the annulus is unstable to density perturbations flLightman fc Eardley 1974 ). This 
region is also thermally unstable (Sunyaev & Shakura 1975). 

As mentioned in §^^, the inner region of a steady thin disc at high M is 
radiation pressure and electron scattering dominated. Hence the region from which 
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X-rays originate is unstable and much effort went into finding how to reconciliate 
this fact with observations. The Shapiro-Lightman-Eardley (SLE) model is one 
such attempt in which the viscously unstable region gives way to a hot, optically 
thin plasma at small radii ( |5hapiro et al. 1976| ). In the SLE model, the protons 
are heated by viscous dissipation and cool by Coulomb interactions with electrons. 
The electrons lose their energy efficiently by local inverse Compton scattering on 
soft photons (explaining the observed hard X-rays) so they have a much lower 
temperature (T c ~ 10 8 K) than the protons (T p ~ 10 11 K). This solution is 
still thermally unstable (Pringle 1976). Lowering the temperature cools the flow 
back to the Shakura-Sunyaev solution. Raising the temperature increases viscous 
heating and decreases the inflow timescale so that protons do not have time to 
give their energy to the electrons before being accreted. Energy is advected and 
the SLE flow becomes a two temperature advection-dominated flow (see §3.2). 

The conclusion that a radiation pressure dominated thin disc is unstable de- 
pends upon the assumptions made. The region can be stabilised by invoking a 
different viscosity law ( |Piran 197§| ), for instance one in which v is proportional to 
the gas pressure only and not the total pressure {fl viscosity). Other possibilities 
include e.g. mass depletion by a wind (Piran 1978) or irradiation heating from the 
boundary layer (Czerny et al. 1986). Most importantly, including radial advection 
of heat stabilizes the flow at higher accretion rates to give a slim accretion disc 
(see §3.2 below and Abramowicz et al. 1988). 



2. 7 The disc instability model 

The stability of a disc can be studied by looking at the local thermal equilibrium 
curves. This is the set of (£, M) such that Q + = Q~ for given M*, a and radius R 
in a steady-state disc. In steady-state, one can use indifferently T c ff or T c instead 
of M (Eq. |2.5|) . If dM / <9£ < somewhere along this curve then the annulus 
is unstable. This does not happen in the Kramers opacity regime (§ [2.6| ). But at 
temperatures lower than 10 K the opacity decreases dramatically due to hydrogen 
recombination. This reverses the slope of the (£, M) relation until the temperature 
has dropped enough that hydrogen is neutral: there is a range of accretion rates 



for which the local annulus is unstable (Meyer & Meyer- Hofmcister 1981). Fig. 1 



shows examples of thermal equilibrium curves calculated at several radii using 
a detailed vertical integration (M* = 7Mq, a = 0.1; the code is described in 



Hameury et al. 1998) 



Such a cu rve, drawing an S in the (£, M) plane, leads to a limit cycle ( Bath & 



Pringle 1982 ). When mass transfer rates M% is in the unstable range, the annulus 
must either jump up on the top branch or down on the bottom branch. On the top 
branch the accretion rate is higher than the average M t so density decreases. When 
the density reaches £ m i n , a slight temperature perturbation causes the annulus to 
cool on a thermal timescale to the lower branch (Q~ > Q + in the region around 
£min)- The accretion rate is lower than M t on this branch so the density increases 
until S max is reached, the annulus heats to the upper branch and the process 
repeats itself. Changes in £ occur on the viscous timescale t v - is oc 1/T C which 
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Fig. 1. Thermal equilibrium curves Q = Q~ calculated at different radii with M* = 
7Mq and a = 0.1. From top to bottom, the S curves are represented in the (S,T C ), 
(E, T ff )m and (£, M) planes (equivalent in steady-state). The middle branch is unstable 
due to hydrogen recombination and the local limit cycle is shown in the middle panel. For 
a large enough disc fuelled by a steady mass transfer rate the inner region is on the hot 
branch while the outer is on the cold branch (easily seen by looking at the intersection 
between the fixed M rmt and the S curves in the last panel). The intermediate region is 
unstable and numerical studies of the non-linear outcome show the disc cycles between 
a hot (outbursting) state and a cold (quiescent) state. 
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is shorter on the cold branch than on the hotter upper branch. The local cycle 
consists of a rapid rise to a brief high M state followed by a longer low M quiescent 
state. In a full disc the situation is complicated since mass is transfered from and 
to each region. The transitions in the (S,M) plane are then more complex than 
the above picture suggests. 

The critical values obtained from vertical structure calculations (Fig. 1) are 
well approximated by fits: 

£ min a 8.3 a-°- 7 M^ A R\i g cm" 2 
S max a 13.4 a-°- 8 Mf - 4 J^6 X g cm" 2 
« 4 • 10 15 a -°Mi°- 9 R% 7 g s" 1 
M Smln a 9 • 10 15 a°- Mr°- 9 Rf 7 g s" 1 

For a steady disc to be everywhere on the hot branch (hence stable) requires 
Mt > Ms mln (^?disc), giving values which are not unrealistic for XRBs and CVs. 
On the other hand, a stable steady cold disc requires M t < -^£ max (-R*) which, due 
to the steep dependence on R, is very low: ~ 10 kg/s for a black hole/neutron 
star and < 10 13 g/s for a white dwarf. Cold, stable steady discs extending down 
to the surface of the compact object are unlikely to exist. 

At intermediate accretion rates, the steady disc solution goes through the un- 
stable branch at some radius. A density front appears a nd th e disc evolves. The 
front speed is set by the radial inflow rate of matter (Eq. 1.11): 

Ur = ~dbA (l/Ei?1/2) 80 «/«^«« c - ( 2 - 14 ) 

since the typical scale of the variations is AR ~ H. The propagation timescale is 
tf = R/vj ~ (ivis^ther) 1 ^ 2 , reflecting that this is a combination of the thermal and 
viscous instabilities. The disc enters a hot state until it cannot sustain the high 
accretion rate needed. A cooling front then propagates, shutting off accretion onto 
the central object and matter piles up until the cycle can start again. 

This disc inst ability model (DIM) has been applied to CVs, XRBs, AGNs and 
protostellar discs ( Hartmann fc Kcnyon 1996 ; Burdcri et al. 1998 ; Lasota 2001 ). S 
curves have been calculated by many different groups using various assumptions 
and degrees of complexity for heat transport or viscous dissipation. Although 
it should really be confirmed from ab initio calculations of viscous transport, the 
reversal of the vT, slope is a robust unavoidable feature. The location of the critical 
E points can vary from model to model: this is equivalent to (small) variations 
in the unknown parameter a. But the critical accretion rates are independent 
of a as (a) H ionisation is what ultimately sets the critical (surface) T c g so it 
is approximately constant (around 8000 K for the hot branch, see Fig. 1); and 



(b) M is related to T e g independently of v (Eq. 2.4). Therefore, the essential 
feature of the DIM is a strong prediction of which systems should be stable. The 
second important aspect is that it offers a framework to understand the outbursts 
of unstable systems. In the following, I discuss the major results obtained from 
the application to cataclysmic variables and X-ray binaries (and defer to Lasota 
200 Qjr a thorough review and complete references). 



14 



Title : will be set by the publisher 



2.8 Models of Cataclysmic Variables 

The DIM reproduces well the distinction between novae-like ( s teady) and dwarf 
novae (unsteady) systems for cataclysmic variables ( Smak 1983| ; |Osaki 19961 ). The 



DIM is als o able to ro ughly reproduce the outbursts of dwarf novae (confirming 



the idea of |Osaki 1974| ) but on one condition: that the viscosity on the hot branch 
is much higher than on the cold branch. Although v varies with temperature, a 
constant a produces only small short outbursts because £ m ; n and S max are too 
close. To explain the amplitudes and timescales, a must change with ahot ~ 



0.1 and a co id ~ 0.01 (Smak 1984). This is one of the most solid conclusions to 
come out of the DIM which can be tested against theoretical models of angular 
momentum transport. In the MM case such a change is not inconceivable as its 
strength depends upon the gas ionisation fraction. Outbursts can therefore give 
some constraints on the viscosity, something steady state models cannot do. 

However, further investigation shows that many additional effects need to be 
considered when modelling outbursts. For instance, the DIM assumes a steady 
mass transfer rate from the secondary but there is no reason for this to be the 
case. In fact, magnetic CVs do show variations of factors 10-100 which can be 
directly attributed to M t since they have no accretion discs. White dwarfs with 
lower magnetic moment can also disrupt the inner regions of a quiescent low M 
disc. Tidal torques from the secondary are important in setting the outer disc 
radius, how it varies and determining the thermal balance in the outer regions. 
The impact of the stream of matter incoming from the Lagrange point and UV 
radiation from the hot white dwarf can also significantly heat the disc and change 



its properties ( Buat-Menard et al. 2001 ). A wide variety of CV lightcurves can 



be accounted for by varying combinations of these effects, using standard values 
for a and some earlier 'failures' of the DIM were actually due to missing some of 



these important phenomena ( Lasota 2001 ) 



However, there are still many unresolved issues and this testifies to large num- 
ber of high-quality observations of CVs. For instance, the DIM cannot readily 
explain the flat quiescent lightcurves or some of the variations observed in the 



outbursts from cycle to cycle (Smak 200C). Moreover, the DIM still fails to make 



accurate spectral predictions, with fits requiring very large values of a. Oversim- 
plification of the radiative transfer may be partly responsible but the issue really 
is that the vertical distribution of the heat generated by angular momentum trans- 
port is unknown. This has a major impact on spectral line formation. A related 
problem not addressed by the DIM is the formation of (line-driven) disc winds 
that can carry away sizeable fractions of the accreted matter. 



2.9 Models of low mass X-ray binaries 

Low mass X-ray binaries are very similar to CVs, with comparable companions, 
orbital periods and Roche lobe overflow mass transfer rates. One would therefore 
expect the DIM to work just as well as in CVs. But the standard DIM wrongly 
predicts that all observed systems should be transient. The solution to this prob- 
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lem is that one cannot neglect irradiation heating in X-ray binaries. High energy 
radiation coming from the inner regions of the flow can interact with the flow 
at larger radii and heat it to the upper branch of the S curve. A small fraction 
(~ 0.1%) of the luminosity deposited this way is enough to make a difference. 
Although there is ample evidence for irradiation in XRBs ( Lewin ct al. 1995] ) , the 
geometry and processes through which high energy photons reach the outer disc 
are still unknown. They could involve a scattering corona, disc warps and/or a jet 
(as in the AGN lamppost model; Dubus et al. 1999). 

A typical dwarf nova has eruptions lasting a few days and recurrence times of 
a few months. A typical transient XRB has an eruption lasting a few months with 
a recurrence period of years ( |Chen et al. 1997a ). Since the main difference with 
CVs is not the mass of the accretor but its radius, one would think that the reason 
for the different outbursts has something to do with the flow reaching in to radii 
10 6-7 cm in XRBs instead of 10 8-9 cm. But numerical simulations show such 
discs undergo large numbers of small amplitude, short duration and recurrence 



time outbursts. This is not surprising as Eq. 2.14 shows that the amount of mass 



needed to trigger an outburst (i.e. for £ > S max ) becomes very small at small R. 
Including irradiation in outburst lengthens the outbursts by keeping the disc hot 
longer but the models are still far from observed lightcurves. Another problem is 
that the accretion rate onto black hole has to be < 10 kg/s in quiescence since 
£ must be lower than S max everywhere in the disc. But black hole XRBs are 
detected at levels of 10 31 erg s _1 which is incompatible with the accretion rate 
and the predicted T e g ss 3000 K in quiescent thin discs. 

As it turns out, there is an elegant solution to both of these problems if the 
thin disc is assumed to be replaced ('evaporated') by a low radiative efficiency 
accretion flow (LRAF) at small radii in quiescence: (1) this removes the highly 
sensitive region which triggered many outbursts; (2) higher accretion rates can 
be 'hidden' in the radiativcly inefficient flow and help increase the recurrence 
timescale (the disc acts as a pierced bucket trying to fill up); (3) spectral models 
of a class of such flows can explain the spectral energy distribution of quiescent 
XRBs. LRAFs are discussed in the next section. 

Models of a modified DIM including irradiation and the evaporation of the 



inner thin disc show good agreement with the observed lightcurves ( Dubus et al 
200lfj~|These models do not require a different viscosity than that of CVs. This is 
gratifying since, after all, the conditions in the disc at large radii are not different 
in CVs and XRBs. On the other hand, the prospect of obtaining strong constrains 
on viscosity by comparing theoretical and observed lightcurves of CVs and XRBs 
fades away as increasingly complex (and poorly known) physics has to be added 
to the model. 



3 Low radiative efficiency accretion flows: thick discs 



The obvious contrast to the thin disc approach is to set / « 1 i.e. to look for 
solutions with little or no radiation cooling to Eqs. 1.9-1.6, The low radiative 
efficiency accretion flow (LRAF) obtained is cooled by the advection of viscous 
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heat (Ichimaru 1977 




Katz 1977 


; iBegelman 1978; 


Rees et al. 1982 




Abramowicz 


et al. : 98S; 


Narayan & Yi 1994 


1995b; 


Abramowicz et al. 1995|). 


Being cooled 



by advection, LRAFs therefore have a radiative timescale which is longer than 
the accretion timescale: t acc /i ra d ~ (1 — f)(Rfl/c s ) 2 <C 1 when / = 1. Inversely, 
this ratio is very large (as expected) for the thin disc (/ = 0) so there is a critical 
parameter value (e.g. M at and a fixed), set by the detailed radiative processes, 
beyond which LRAFs can exist (see §3.2). 

LRAFs are necessarily hot flows with T < T g (see Eqs. 1.18-1.19). A special 
case is that in which / = 1 and fl — which corresponds to a purely radial 
accretion flow i.e. to the Bondi-Hoyle-Lyttleton solution describing accretion from 



the interstellar medium (Bondi & Hoyle 1944). In this sense, the adiabatic flows 
discussed below are a generalisation of this solution to the case of non-zero initial 
angular momentum. 

Because of their high temperatures and low radiative efficiencies LRAFs pro- 
vide an ideal context to model low luminousity objects showing emission at high en- 
ergies (e.g. Lasota et al. 1996 ) and have thus attracted considerable interest. Since 
T < T g implies H < R, adiabatic flows are geometrically thick and more amenable 
to numerical simulations than thin discs in which radiation processes have to be 
taken into account and where resolving the very small vertical scale height is nu- 
merically challenging. Theoretical work on LRAFs has thus been enriched by the 
possibility to simulate adiabatic accretion, including ab initio transport of angular 
momentum via the MRI (hence requiring no additional assumptions) . 

The equations for a steady state accretion flow with an a type viscosity can be 
put in dimensionlcss form as shown in §1.6. Far from the boundary conditions, the 
equations have no length scale and a self-similar solution in which the variables are 
written in powers of R is possible. The thin disc is one special case with / < 1 and 



the ne xt se ction describes the extension to arbitrary / (spruit et al. 1987; Narayan 
& Yi 1 994 ). In principle / is found self-consistently from the radiation transfer 
and this is done in a rudimentary way in §3.2, leading to the "unified description of 
accretion flows". Some applications are presented in §3.3 while current problems 
and controversies in LRAFs are summarized in §3.4. 



3.1 The self-similar solution 



Assuming / is constant, Eqs. 1.18- 1.19] immediately show that O oc Six oc i?~ 3 / 2 
and T cx T g oc i2 _1 . Using these in Eq. 1.17 gives v r oc Rr 1 / 2 . The hydrostatic 
balance HQk — c s oc T 1 / 2 gives H oc R, mass con servation gives S oc R~ x l 2 so 
p oc R~ 3 ' 2 . With these dependences, Eq. 1.18 - 1. IS solve for ft and T from which 
the other variables can be deduced: 



n 2 

v r 
T 

H/R 



(e/f) h(a n j) n 2 K 
-(3/2)a h(«»7,/) RM K 
(3/2)(l + e) h(a )7 J) T s 

/A 



(3.1) 
(3.2) 
(3.3) 
(3.4) 
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where e = (5/3 — 7)/(7 — 1) and h is a function^ of a, 7 and /. To fix ideas, 
fi w 0.7S1 K , Vr ~ 5- 10 7 cm-s" 1 , T w 4- 10 9 K and ff/i? w 0.75 at 10 9 cm around a 
10 M black hole with / = 0.75, a = 0.1 and 7 = 4/3. The flow is sub-Keplerian, 
very hot and the accretion timescale is very short compared to the thin disc. 

The thin disc is recovered as expected when / — > 0, giving h r* f/e — ► 0. In 
the other limit, as / increases towards 1, h increases so the temperature, viscosity 
and inflow velocity v r ~ i//i? all increase. Since i acc = R/v r ~ (l//i)ithor, the 
radial inflow time can become comparable to the thermal time. The disc becomes 
sub-Keplerian and thick with H/R < 1. If the gas is adiabatic with 7 = 5/3, the 
flow becomes purely radial with O — ^ i.e. the solution tends to Bondi spherical 
accretion. Ogilvie (1999| ) studied the evolution of an initial torus of gas with some 
angular momentum at R accreting adiabatically, a setup similar to that described 
for thin discs in §2.5. He finds the self-similar solution with / = 1 is asymptotically 
approached with time. 

There are some obvious limitations to the self-similar solutions. The self-similar 
solution extends at all R and does not address what happens at boundaries. But 
integration of the equations with physically motivated boundary conditions show 
the self-similar solution describes the flow well enough far from the boundaries 
( |Chen et al. 1997b ; Narayan ct al. 1997). Global solutions to the equations in 
gene ral relativity have also bee n found ( Abramowicz et al. 1996 ; Pcitz fc App] 
1997; Popham & Gammie 1998). The validity of using the vertically integrated 
set of equations is not obvious when H/R ~ 1. [Narayan fe Yi (1995a| ) studied 
self-similar solutions to the equivalent 2D (R,8) set of equations and found that 
the above solution was a good approximation if the integration is interpreted in 
terms of 9 instead of z. 



3.2 Radiative processes: ADAFs and slim discs 

The parameter / was assumed given when it should really be found self-consistcntly 
from the radiative processes at all radii. In practice, this implies guessing /, 
determining Q~ for all relevant radiative processes (e.g. optically thick flux or 
free free, synchrotron, pair creation, comptonisation etc.) using the values for E 
and T from the radial solution, work out the improved value for / and repeat the 
process until convergence is reached (e.g. Narayan fc Yi 1995b] ). To complicate 



things further there is no a priori reason for / to be constant in the flow but this 
is an assumption which is often made for simplicity and is likely to be reasonable 
in light of other approximations in the radiative transfer ( Esin et al. 1997] ). 



3.2.1 Optically thin case: ADAFs 

Here, all of these complications are left aside and simple rudimentary forms for 
Q~ are used to find self-consistent solutions. The energy equation (Eq. [1.12 ) in 



1 h(a, 7, /) is the solution to 9a 2 h 2 + 4(5 + 2j)h - 8 = 
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steady state can be rewritten as (Abramowicz ct al. 1995) 




Qadv = Q+ is - Q" - * (^j = [ £ ] + U Q- (3.5) 

where u = dlnO/dmii and £ = dbx pT 1 / 1 -^ /din R + 4(1 - /?)dlnpT- 3 / dln -R 
for a mixture of perfect gas and radiation (/3 is the ratio of gas pressure to total 
pressure). In an adiabatic flow Q~ = so the equation reduces to £ = (j(RQ/ 'c s ) 2 . 
Adiabatic accretion is possible for any M when this is verified. 

If Q~ is not negligible and the flow is emitting optically thin free-free ra- 
diation then Q~ ~ HplT 1 / 2 ~ tt K T, 2 (J3 = 1). For a given R, £ and a the 



mass accretion rate M satisfying the thermal equation (Eq. 3.5) depends on 



global properties of the flow. The thermal equation can be made local by as- 



suming that the self-similar scalings hold and that ft — Qk- Eq- 3.5 then be- 
comes M 2 - adERi/^M + aC 2 S 3 i? 2 = (d, C 2 constants). Below a criti- 
cal E cr ; t cx a/R there are two possible solutions for M: one with Q~ > Q 



adv 



corresponding to a Shapiro-Lightman-Eardley flow (see §2.6) and the other with 
Qadv > Q is called the advection-dominated accretion flow (ADAF) in the lit- 
erature. Above the critical M cr it oc a 2 i? -1 / 2 corresponding to S C rit there is no 
solution. M cr it can also be found by equating the accretion timescale to the radia- 
tive timescale for bremstrahlung: ADAFs appear when the latter becomes longer 
than the former. 

3.2.2 Optically thick case: slim discs 

In the optically thick case we have already seen that one solution is the Shakura- 
Sunyaev thin disc solution (Q + = Q~). As M increases, radiation pressure in- 
creases and the solution becomes unstable (§|2.6]). Writing the cooling term in the 
Eddington approximation with electron scattering opacity and assuming radiation 
pressure dominates, the thermal equilibrium can be rewritten with H/R instead 
of M as: 

H\ 3 fH\ Q 2 C 



Qad V = Q v + is -Q-^j =[x)^-^M^ (3-6) 

where C is a constant. Again, this equation can be made local by assuming 
Keplerian rotation and £ constant. For a given E, R and a there can be up to 
three solutions for H/R. Detailed calculations show an S-curve in the (Y,,H/R) 
plane where the bottom branch is the Shakura-Sunyaev solution, the middle branch 
the unstable P rac i dominated thin disc and the top branch a new solution called 
the slim disc solution. Slim discs are radiation-pressure dominated and cooled by 
advection. 

3.2.3 Thermal equilibrium curves 

Both the ADAF and slim disc branches are thermally and viscously stable (see 



Honma et al. 1991; Kato et al. 1997, and references therein for details). An easy 
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Fig. 2. A global view of the thermal equilibrium curves in the column density E, mass 
accretion rate rh (normalised to the Eddington rate) plane. This figure, adapted from 
Chen et al. (1995), shows the branches corresponding to different types of accretion. 
The curves are calculated at a radius of 30 R g for a black hole of 10 M© and using 
a — 0.1 (green/light lines) or a = 1 (red/dark lines). The dark areas enclose the regions 
where Q + > Q~ in the case a = 0.1 and Q~ > Q + in the case a = 1. This shows 
the SLE solution (dotted lines) is thermally unstable in both cases (a slight increase in 
temperature at constant E pushes a disc on the SLE branch in the Q + > Q~ region so 
that it cannot cool). On the other hand, the ADAF, slim disc and gas pressure dominated 
thin disc branches are stable. The radiation pressure dominated thin disc (dotted line) 
is thermally and viscously unstable. 



way to see this is to plot the different branches found above in the (E,M) in a 
similar fashion to what was done for the S-cur ves of the disc ins tability model in 
§|/| This is shown in Fig. 2 (adapted from |Chen et al. 1995| ). At low E, the 



disc is optically thin and there are two solutions: the thermally unstable SLE 
solution and the stable ADAF branch. At higher E, the disc is optically thick and 
there are three branches: the Shakura-Sunyaev thin disc (i.e. the upper branch 
of the S-curve in Fig. 1; the S corresponding to hydrogen ionisation occurs at 
much lower values of M than shown in Fig. 2; see §2.7.), the radiation pressure 
dominated thermally/ viscously unstable thin disc and the slim disc branch at 
high M . There is a change of topology at a critical a above which an advection- 
dominated solution is possible for all M. Although these curves are a useful tool 
to illustrate the various accretion solutions, it is worth stressing that they are local 
representations of an inherently non-local phenomenon (advection) and that they 
assume an a type viscosity. These curves therefore are not in any way a complete 
description of accretion flows. 



3.3 A p plica tions of L RA Fs 

LRAFs have applications to a variety of different objects in various situations. The 
S-curve at the transition between a thin disc and a slim disc (see Fig. 2), when 
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radiation pressure takes over, opens up the possibility of a thermal-viscous cycle 



like that described in §2.7. The associated mass accretion rates are very high and 
the timescales are very short so the cycle timescale is fast, of the order of seconds to 
minutes. This has found possible applications in some XRBs which accrete at high 



rates (e.g. Cannizzo 1996; 3zuszkiewicz & Miller 2001, and references therein) 



Optically thin LRAFs have many more applications. As mentioned at the 
end of §2, models of transient XRBs suggest a LRAF replaces the thin disc at 
small radii in quiescence. Support for this comes from models of optically thin 
LRAFs that reproduce the spectral energy distribution of quiescent X-ray binaries 



(Narayan et al. 1996). A comparison of the quiescent luminosities of neutron 
star transient XRBs against black hole transient XRBs shows that the neutron 
star systems are systematically more luminous by an order of magnitude or two 
than their black hole counterparts. This has been interpreted as evidence for a 
horizon in black holes. In both cases the compact object is surrounded by a LRAF. 
However, the energy stored in the flow is radiated when it reaches the surface of 
the neutron star while it disappears through the hori zon in a black hole, explaining 



the luminosity difference (for a recent discussion see Garcia et al. 2001 ) 



The Esin et al. (1997| ) model of transient XRBs spectral states is a natural 



extension of the idea that LRAFs and thin discs coexist. As the transient XRB goes 
into outburst, the ADAF retreats to smaller radii and eventually is fully replaced 
by the thin disc when the increasing mass flow rate is above the ADAF critical mass 
accretion rate (§3.2). Since the thin disc has a much lower temperature than the 
ADAF this roughly implies a transition between a hard and a soft X-ray spectrum 
during the rise to outburst and vice-versa during the return to quiescence. These 



ideas have been applied to the observed spectral changes in Nova Mus 1991 (Esin 



et al. 1997) 



However, the most studied application of LRAFs is to Sgr A*, the black hole at 
the center of our Galaxy. Despite its dynamically- measured mass of a few 10 6 M , 
the X-ray luminosity of Sgr A* is only 10 33 ergs s _1 , well below the Eddington 
luminosity. The mass accretion rate inferred from observations is much higher than 
the luminosity suggests making it an ideal candidate for a LRAF. The dormant 
massive black holes at the center of n earby galaxies are also prime candidates for 
LRAFs (for a review of the above see Narayan 2002] ) . 



3.4 Current issues 

3.4.1 LRAFs around neutron stars and white dwarfs 

Matter accreting onto a black hole passes necessarily through the last marginally 
stable circular orbit and a sonic point before free-falling through the horizon. The 
inner boundary condition of global LRAF solutions around black holes is therefore 
provided by the regularity requirement at the sonic point. This sets the specific 
angular momentum accreted by the black hole. In contrast, matter falling onto 
a neutron star will have a vanishing radial velocity and may not necessarily pass 
through a sonic point (this depends on the neutron star radius i.e. equation of 
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state). The subsonic flow can spin-up or spin-down the neutron star. Energy 
released at the surface of the neutron star can significantly change the luminosity, 
spectra and heat budget of the accretion flow. All of these features make accretion 
onto a neutron star (or white dwarf) a difficult problem to tackle (see Medvedev 
& Narayan 2001, and references therein). 



3.4.2 Two temperature plasmas 



Optically thin LRAF solutions exist only below the critical mass accretion rate 
at which the accretion timescale is shorter than the radiation timescale (§3.2). 
Since it is mostly the electrons which radiate, this mass accretion rate is set by 
the electron temperature. To be astrophysically interesting, the electrons must be 
assumed to have a much lower temperature than the protons (Shapiro et al. 1976; 
Rccs ct al. 1982|) . This is achieved if viscous heating mostly goes to the protons 



and if these exchange heat with electrons only via Coulomb collisions (another is- 
sue is whethe r protons and electrons have e nough time to thermalize before being 
accreted, see Mahadevan fc Quataert 1997 ). All applications have assumed such 
two temperature LRAFs (the self-similar solution of §3.1 assumes a single tem- 
perature for simplicity). The relative fraction S of heating t hat ultimately goes to 
the electrons depends on the microphysics and is uncertain (Narayan & Yi 1995b; 
Quataert 1998| ; |Quataert fc Gruzinov 1999 ). Plasma i nstabilities could erase any 



temperature difference between electrons and protons (Bcgclman & Chiueh 1988 
Bisnovatyi-Kogan fc Lovelace 1997 ) . Observational constraints can in principle be 
obtained since the radiation spectrum depends on 6. 



3.4.3 ADIOS 

Another interesting feature appears when looking at the total local energy of the 
gas (per unit mass), defined as the sum of the kinetic, internal and potential energy 
from gravitation and pressure: 

1, , 99n GM P 3e(/- 1/3) , , , , 
E = e c + e l + e g + e p = -{v 2 r +R 2 n 2 )+e- — + -= ^ - ± > R 2 tf K (3.7) 

This is also the Bernouilli integral, which is constant along streamlines for an 
adiabatic inviscid fluid (a=0 and / = 1). Since the number is positive for / > 1/3, 
a particle going outwards on a streamline could reach infinity with a net positive 
energy, i.e. the fluid is energetically unbound and there could be an outflow. The 
local energy is positive because potential energy liberated at small radii is radiated 
at larger radii just like in the thin disc. When the flow cools efficiently this extra 
energy is efficiently radiated locally and the Bernouilli number is negative (which 
explains the / = 1/3 limit). Note that in the self-similar solution there is an 
infinite amount of energy released at R = to redistribute hence it is unclear 
whether realistic boundary conditions will also yield a positive Bernouilli integral 
(Abramowicz et al. 2000) . 

The assumption that radiatively inefficient flows will necessarily lead to out- 
flows is referred to as the Advection Dominated Inflow-Outflow Solution (ADIOS) 
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after Blandford fc Begclman (1999 ). Work on ADIOS has focused on extensions 
of the ADAF with an outflow parameterised by M oc R p (p constant). Whether 
Nature prefers ADIOS to ADAFs is an open issue with no clear observational tests 
( |Quataert fc Narayan 1999 ) and conflicting numerical simulations ( |Hawley et al 
200 1|; feumenshchev fc Narayan 2002[). 



3.4.4 CDAF 



It may turn out that neither ADAFs or ADIOS are correct descriptions. ADAFs 
are unstable according to the Hoiland criterion for convective stability of a rotating 
flow: 



stable if (1 



J- a 



dlnP\ dln(P 1 /T/p) 



dln.fi 



dlnfi 



? dln(fi 4 f7 2 ) , . 



where the first term is the (squared) Brunt- Vaiasala frequency (the frequency at 
which a perturbed blob of gas oscillates around its position in a stratified atmo- 
sphere) and the second term is the (squared) epicyclic frequency (the frequency at 
which a perturbed particle oscillates around its orbit). Convection will transport 
angular momentum and energy, changing the radial structure of the flow. The 
Convection Dominated Accretion Flow (CDAF) is one solution which assumes 
convection will transport angular momentum inwards (by analogy with axisymet- 
ric hydrodynamic flows) and that the steady state flow is close to the marginally 
stable limit at which convective angular momentum transport balances outwards 
viscous transport. One condition for this to happen is that viscous transport is 
not too strong, i.e. for low values of a. There is then a net flux of energy out- 
wards but very little net accretion (as opposed to ADAFs), yielding the CDAF 
self-similar scaling p oc Rr 1 / 2 (Narayan et al. 2000; Quataert & Gruzinov 2000). 



CDAFs have attracted a lot of attention because this density scaling was found in 
several numerical simulations (e.g. Stone et al. 1999). CDAFs are highly debated 
for several reasons: (1) purely hydrodynamical simulations with a Navier-Stokes 
viscosity may not catch the basic behaviour of a realistic flow; (2) MHD simu- 
lations may be inappropriate because the plasma is collisionlcss; (3) the Hoiland 
criterion may be inappropriate; (4) CDAFs may be thermodynamically impossible 
since energy dissipated in the MRI turbulent cascade cannot be used to generate 



convection (see e.g. Balbus & Hawley 2002; Hawley & Balbus 2002; Igumenshchev 



fc Narayan 2002; Narayan et al. 2002) 



4 The transition between thin and thick discs 

As illustrated by Fig. 2, accretion with a given M has a choice between several 
different types of flows. There is at present no guidance as to which particular 
type of flow is preferred, which admittedly an important open question when 
considering the potential applications of LRAFs. The flow may always prefer 
a LRAF whenever it is possible and this is called the strong LRAF (or ADAF) 
principle; or it may choose the LRAF only when there is no other solution available 
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(the weak principle). One application in which this has important consequences 
is the LRAF+thin disc model of transient XRBs. Insights into this issue might 
be gained by studying the directly related problem of the physics of the transition 
from a thin disc to a LRAF which is briefly discussed here. 

4.1 Condition for the transition 

Despite its shortcomings, an examination of Fig. 2 suffices to show most of the 
problems linked with the transition from a LRAF to a thin disc. A first possibility 
is to have a slim disc - thin disc transition but this occurs at high M and there will 
be an unstable region (see §3.3). When a > a cr ;t the thin disc makes a transition 
to the Shapiro-Lightman-Eardley type flow rather than the slim disc. However, 



the SLE flow is unstable (£2.6). These are the only possible smooth transitions 
with a thin disc and both will are unstable. The transition from a thin disc to 
an optically thin LRAF (ADAF), which is observationally motivated, necessarily 
involves some discontinuity in the framework of Fig. 2. At the transition, flow 
variables such as the radial inflow speed v r or the temperature T must increase 
by several orders of magnitude and this involves some additional input of energy 
which is not taken into account in Fig. 2. 

4.2 Possibilities 

Studies of the transition between an optically thin LRAF and a thin disc therefore 
always involve some additional energy source. In the vertically averaged model 
of Honma (1996| ) the extra term comes from turbulent energy transport. This 



term dominates in the transition region. The transition is likely to be a 2D or 3D 
process that occurs gradually in radius i.e. the thin disc evaporates into a more 
tenuous gas in a sandwich around it. The study is therefore not much different 
from the studies of disc corona models. The generic idea is that radiative cooling 
is oc p 2 so that any heating of the tenuous upper layers of a thin disc can lead to 



runaway heating (see §2.6) and from there to conditions reminiscent of LRAFs. 



The energy terms that can be involved include viscous heating (usually taken to 



be oc P hence oc p, §|2.6| ), heating by the reconnection of magnetic loops (Galeev 
et al. 197Sj ) , irradiation by photons or ions from the inner flow (Deufel & Spruit 



200C), or electron conduction ( Meyer et al. 2000| , and references therein). Progress 



requires a better understanding of how the energy is dissipated, transported and 
radiated, much of which is based on difficult microphysics. This issue has started 
to be tackled by MHD numerical simulations ( Miller fc Stone 200C| ). 
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